How Gravitational-wave Observations Can Shape 
the Gamma-ray Burst Paradigm 



I. Bartos 1 , P. Brady 2 , S. Marka 1 

1 Department of Physics & Columbia Astrophysics Laboratory, Columbia University, 
New York, NY 10027, USA 

2 Center for Gravitation and Cosmology, University of Wisconsin-Milwaukee, 
Milwaukee, WI 53211, USA 

Abstract. By reaching through shrouding blastwaves, efficiently discovering off-axis 
events, and probing the central engine at work, gravitational wave (GW) observations 
will soon revolutionize the study of gamma-ray bursts. Already, analyses of GW 
data targeting gamma-ray bursts have helped constrain the central engines of selected 
events. Advanced GW detectors with significantly improved sensitivities are under 
construction. After outlining the GW emission mechanisms from gamma-ray burst 
progenitors (binary coalescences, stellar core collapses, magnetars, and others) that 
may be detectable with advanced detectors, we review how GWs will improve our 
understanding of gamma-ray burst central engines, their astrophysical formation 
channels, and the prospects and methods for different search strategies. We place 
special emphasis on multimessenger searches. To achieve the most scientific benefit, 
GW, electromagnetic, and neutrino observations should be combined to provide greater 
discriminating power and science reach. 
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1. Introduction 

Gamma-ray bursts (GRBs) are the brightest electromagnetic explosions in the Universe 
(e.g., [1]). For the short time of their activity, they outshine all other sources in the sky 
in gamma rays. These energetic explosions originate from cataclysmic cosmic events, 
whose "inner engine" that drives the observed emission is confined to volumes mere 
tens of kilometers across, as indicated by the duration and variability of gamma-ray 
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emission (e.g., [2]). This inner engine is hidden from direct electromagnetic observations, 
even though the variability of the emission is directly related to the engine's activity. 
Gravitational waves (GWs), on the other hand, carry information directly from the 
inner engine to the observer. 

With the near completion of advanced GW detectors [3, 4, 5, 6], we will soon be 
able to directly probe the dynamical processes leading to the creation of GRBs, which 
are less (or not) accessible via other messengers. The Advanced LIGO [3] and Virgo 
[7] detectors are planned to begin observation in 2015, albeit initially below design 
sensitivity, gradually reaching their design sensitivity around the end of the decade 
[8]. Upon reaching their design sensitivity, the Advanced LIGO and Virgo detectors 
are expected to be about 10 times more sensitive (at their most sensitive frequency 
band around ~ 150 Hz) than the initial LIGO/ Virgo detectors. Further, advanced 
detectors will have a wider sensitive frequency band, reaching down to ~ 10 Hz. This 
will improve their sensitivity to wide-band sources such as compact binaries. The global 
GW detector network will further increase its reach in both distance and source direction 
reconstruction by the addition of further detectors [9], such as the Japanese KAGRA 
(formerly LCGT; [6]) detector that is under construction, as well as the planned third 
LIGO observatory in India [8] . 

There are a number of mechanisms conjectured to result in a GRB. Potential 
mechanisms include the formation of a central object surrounded by an accretion disk 
from the merger of binary neutron stars [10, 11, 12] or a neutron star and a low-mass 
black hole [13], the collapse of the core of a massive star [14, 15, 16], the global 
reconfiguration of the magnetic fields in magnetized neutron stars [17, 18, 13], the 
accretion-induced collapse of white dwarfs following accretion from a non-degenerate 
companion [19, 20] or following the merger of a white dwarf binary [21, 20], or even 
cosmic strings [22] . Identifying the mechanisms behind the creation of individual GRBs 
will be greatly aided by GW detections. 

This review intends to survey the prospects of GW measurements in understanding 
the central engines of GRBs. While our theoretical view of GW emission mechanisms 
connected to GRBs is rapidly evolving, the main directions, the opportunities and 
limitations of observational GW astrophysics are becoming clearer. With this review we 
aim to provide guidance to (i) astronomers interested in utilizing the capabilities of GW 
measurements to complement electromagnetic or neutrino observations in understanding 
GRBs, and (ii) GW scientists interested in the open questions in GRB astrophysics that 
can be addressed through GW measurements. 

Compact binary mergers, the likely progenitors of most short GRBs (e.g., [13, 23], 
although see [19] for alternative models), are one of the primary targets of GW searches 
(e.g., [24]). Advanced GW detectors will be able to, among others, measure the 
properties of the binary [25], directly probe the dynamics of the central engine, and 
correlate the binary's masses and orientation with prompt and delayed electromagnetic 
counterparts. These observations will constrain the central engine mechanism, provide 
standard sirens for measuring cosmological parameters [26, 27, 28], and constrain the 
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nature of nuclear matter (see Section 3.1). 

GW emission from isolated central engines, such as stellar core collapses, accretion 
induced collapses or magnet ars, is typically at higher frequencies and at lower amplitudes 
than for compact binaries. GW emission depends on the development of varying 
quadrupole moment, which is highly source and model dependent. The detection of 
GWs from isolated objects would help us understand the internal processes in massive 
stellar collapses, including the development of differentially rotating protoneutron stars, 
the nature of nuclear matter, or the creation of massive accretion disks (Sections 3.2, 
3.3 & 2.4). 

For many topics touched upon here, we can refer the interested reader to existing, 
excellent reviews; on GRB blastwaves [1, 13, 29, 30, 31]; on short and long GRB 
observations [13, 32]; on models of central engines involving core-collapse events [33], 
neutron star mergers [34, 35], and isolated neutron stars [36, 37]; on GW emission 
from core-collapse events [33, 38] and binary mergers [39, 40, 41, 42, 43]; on short 
GRB rates [13], long GRB rates [44, 45], and compact object merger rates [46]. 
Some investigations of isotropic prompt and delayed electromagnetic counterparts have 
recently been summarized in [47] (for short GRBs). 

This review is organized as follows. Section 2 outlines GW emission processes 
from GRB progenitors, focusing on the emission scenarios that may produce signals 
which could be detectable with advanced GW detectors. Section 3 presents some of the 
intriguing astrophysical questions that could be answered through the detection of GWs 
from GRB progenitors, in some cases in coincidence with other messengers. In Section 
4 we review GW observation strategies from single GW to multimessenger searches, 
also outlining prospects for the advanced detector era. Finally, Section 5 we briefly 
summarizes the presented results and the near future of GRB astrophysics with GWs. 

2. Gravitational- wave emission processes in GRB progenitors 

In this section, we outline the main GRB progenitor models with an emphasis on their 
expected GW signature. It should be clear that, depending on the progenitor model, 
the GW emission from GRB progenitors varies widely with respect to predicted signal 
strength and characteristic frequency. 

2.1. Compact binary coalescence 

The majority of short-hard GRBs are thought to be powered by the merger of NS-NS 
or NS-BH binaries [48, 11, 12, 34, 42, 49, 50, 51]. Compact binaries are also strong 
GW sources: the GWs produced in the minutes before the merger should be detectable, 
by second generation instruments (i.e., Advanced LIGO and Virgo), reliably out to 
~ 450 Mpc for NS-NS binaries, and further for NS-BH systems [46, 52]. 

Theoretical studies provide detailed information about the binary evolution and the 
gravitational waveform expected from these systems. Moreover, numerical simulations, 
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which include the physics of the NS matter, continue to improve (see Refs. [35, 41, 42, 43] 
for reviews). Here we present an overview of compact binary evolution, which is 
summarised in Fig. 1. At the highest level, there can be four phases: the inspiral 
phase during which the orbit shrinks due to energy and angular momentum loss via 
GWs; the merger phase during which dynamical effects and microphysics of the NS 
become important; the accretion phase when material stripped from the NS(s) accretes 
onto the merged objects; and the ringdown phase when the merged objects settle down 
to an unperturbed (Kerr) BH or a NS. The details vary by NS-NS and NS-BH merger 
and by mass (and other parameters) of the merging objects. The accretion stage is 
likely essential to the formation of short GRBs from binary coalescences, although 
the formation and evolution of protomagnetars following NS-NS merger has also been 
proposed as a plausible central engine [19]. Figure 2 presents the main features of GW 
emission in a schematic spectrum of the GW effective amplitude, separately for NS-NS 
and BH-NS mergers. Below we discuss the different stages of the evolution in more 
detail. 

2.1.1. Inspiral phase — Early in the evolution of a compact binary system, the two 
objects are separated by a relatively large distance compared to their radii. The 
binary elements spiral towards each other by losing angular momentum via the emission 
of GWs. During the early stages of the inspiral, the two compact objects can be 
approximated as point masses for the purposes of GW emission [43] . At the late stages 
of the inspiral, the internal structure of the objects become increasingly important. 
For example, the tidal deformation of NSs in a binary system can (slightly) affect the 
orbital period (and therefore the gravitational waveform) in the late inspiral phase 
[53, 54, 55, 56, 57]. Further, general relativistic spin-spin or spin-orbit coupling can 
cause the binary's orbital plane to precess, affecting the binary's evolution and GW 
emission [58, 59, 60]. 

Nevertheless, the dominant features of the GW signal from the inspiral phase are 
captured by neglecting the spins and internal structure of the binary elements. As the 
objects spiral together, their orbital frequency increases producing a GW signal that 
sweeps upward in frequency. About ~ 15 minutes before merger, the GW from the 
inspiral of a NS-NS binary begins to sweep upward from ~ 10 Hz through the band 
of Earth-based GW interferometers. The effective amplitude h e s = f\h(f)\ of the GW 
signal from a binary system decreases as h e s oc f~ x ^ [61], up to a mass-dependent cutoff 
frequency / cut ~ 1 — 3 kHz [62, 63, 61]. The frequency ranges < 1 kHz and 1 — 3 kHz are 
traditionally considered the inspiral and early merger phases, respectively. For / < f cut 
the merger retains a binary-like structure and consequently emits relatively strong GWs 
[61]. 

Advanced detectors will be able to detect a NS-NS inspiral up to ~ 450 Mpc, 
while NS-BH inspirals will be detectable up to ~ 950 Mpc [46] (the distances are 
given for untriggered searches, with optimal source orientation and direction; for further 
details see Section 4.1.1). The effective survey volume determined by averaging over sky 




Figure 1. Schematic diagram of the evolution of compact binary coalescences. The 
frequency of the emitted GW is indicated for the different stages. NS-NS inspirals are 
observable for a few seconds to minutes. Upon the merger of the NSs, a binary with 
total mass Mu nary > 3 M promptly collapses into a BH. For non-equal-mass binaries, 
the forming BH will be surrounded by an accretion disk. NS-NS binaries with total 
mass MNS,max < ^binary < 3 M (where MNS,max is the mass limit of non-rotating 
NSs) form a hypermassive NS with strong differential rotation, which assumes a non- 
axisymmetric ellipsoid shape. The hypermassive NS survives for milliseconds to a 
second, eventually collapsing into a BH, potentially with an accretion disk. Very low 
mass NS-NS binaries (Mbi nary < Mjsrs,max) can leave a stable NS behind. For BH-NS 
binaries, after an inspiral phase observable for seconds to minutes, the NS either gets 
tidally disrupted (if tidal disruption at radius Rudai occurs before the NS could reach 
the ISCO at Risco), or it plunges into the BH (if Rudai < Risco)- Tidal disruption 
results in a BH with an accretion disk, while no accretion disk forms upon plunge. This 
merger phase, along with the ringdown of the BH after plunge, lasts for milliseconds. 

location and inclination of the sources is ^ 4^(^/2. 26) 3 /3 [46]. Using the current best- 
guess rates of mergers, this gives tens of NS-NS and a few NS-BH binaries detected with 
advanced detectors each year [46]. Additional advanced detectors, such as KAGRA [6] 
or LIGO India [64] , can significantly increase this range [9] . Third generation detectors 
are expected to reach an order of magnitude farther than advanced detectors, i.e. to 
several Gpc, and hence will be able to observe tens of thousands of events a year (e.g., 
[65]). 

2.1.2. Merger phase — Depending on the binary system, the merger can progress in 
multiple distinct directions with qualitatively different GW and gamma-ray emission. 
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Figure 2. Schematic spectrum of effective amplitude h e & for compact binary 
coalescences, (left) NS-NS binaries. During the inspiral phase up to ~ 1 kHz and the 
early- merger phase up to / cut ~ 3 kHz, the system retains its binary-like structure and 
/i e ff scales as / _1 / 6 . If a BH is promptly formed, matter quickly falls in the BH, losing 
angular momentum through emitting GWs around a peak frequency / pea k ~ 5 — 6 kHz. 
If a protoneutron star is formed from a NS-NS binary, it will radiate GWs through its 
quasiperiodic rotation at / qp d ~ 2 — 4 kHz. After matter falls into the BH, the BH rings 
down, emitting GWs at « 6.5 — 7 kHz with exponentially decaying amplitude, (right) 
BH-NS binaries. During the inspiral phase, /i e ff scales as / _1 ^ 6 . If the NS is tidally 
disrupted before reaching the ISCO, GW emission will cut off at /tidal ~ 2 — 4 kHz, 
i.e. the GW frequency at tidal disruption [66, 67]. If the NS plunges into the BH 
without being tidally disrupted, the plunge cuts off GW emission from the binary and 
excites the quasinormal mode of the remaining BH, which rings down emitting GWs 
at frequency / qn m (NS-NS representation was partially inspired by Kiuchi et al. [61]; 
BH-NS representation is based on Kyutoku et al. [66]). 

The formation of a massive accretion disk is probably crucial to the generation of a 
GRB from binary mergers^ [42] . Accretion disks can be created via the tidal disruption 
of a NS at some point during the merger [72, 73, 74, 75, 76, 77]. Alternatively, for NS- 
NS mergers, material with centrifugal support can be left outside the newly formed BH. 
Whether an accretion disk is formed or not depends on the binary properties (mass, spin, 
etc.), as well as the NS equation of state (EOS). A disk mass of ~ 0.01 M , where M is 
the mass of the sun, is probably sufficient to supply the energy for the creation of a short 
GRB [61, 51]. In a recent comparison between numerical simulations and observation, 
Giacomazzo et al. [51] found that the observed emission of short GRBs implies torus 
masses < 0.01 M , which favors "high-mass" NS-NS mergers (M tata i > 3M ). Further, 
they find that BH-NS mergers, while cannot be excluded, would require a very rapidly 
spinning BH (with spin > 0.9). 

Here we list the possible outcomes of the merger and outline the scenarios that can 
lead to them. 

X Millisecond magnetars represent plausible alternative central engine candidates for some GRBs 
[68, 19, 69, 70]. See also, e.g., [71] and references therein, for further alternative models. 
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• BH with no accretion disk — for BH-NS binaries, if the radius of the innermost 
stable circular orbit (ISCO) is greater than the tidal disruption radius, the NS plunges 
into the BH before it could be tidally disrupted, resulting in a BH with no accretion 
disk. This will be the case for binaries with relatively high BH:NS mass ratio (> 5 : 1) 
[74, 75, 42, 78]. This ratio strongly depends on the BH spin (due to the change in the 
location of the ISCO) and the NS EOS [76, 66, 79, 80]. 

For NS-NS binaries with equal masses, if the binary mass exceeds a threshold M t h rj 
the NSs will promptly collapse to a BH upon merger [81, 82], leaving essentially no 
accretion disk behind [73, 83, 82]. M t h r depends on the NS EOS. Kiuchi et al. [55] 
found that a BH is promptly formed if total mass of the binary system is > 3M . 
In this scenario no accretion-powered GRB will be produced, although there are 
possible channels through which even such a system can emit electromagnetic 
radiation [84, 85]. 

• BH with accretion disk — For BH-NS binaries with relatively low mass ratio 
(< 4 : 1), the NS will be tidally disrupted before falling into the BH, which leads to 
the formation of a massive accretion disk [86, 74, 75, 42]. For spinning BHs, the mass 
ratio below which disruption occurs is even higher [76, 78, 87, 66], and the mass of the 
formed disk can greatly depend on BH spin and spin-alignment [76, 88, 66, 80]. Prior 
to tidal disruption, the emitted GW frequency reaches /tidal ~ 2 — 4kHz§ [66, 67]. 
For NS-NS binaries with unequal NS masses and sufficiently large total mass (> 3M 
[55, 83, 90]), the less massive NS will be tidally disrupted, followed by the more massive 
NS's prompt collapse into a BH, leaving a potentially massive accretion disk behind 
[81]- 

In this scenario, after the merger of a NS-NS binary, if the forming disk features 
azimuthal variations, GWs may be emitted by the material orbiting the central object 
around a peak amplitude / pea k ^ 5 — 6 kHz [61]. For BH-NS binaries, / pea k can be 
significantly lower, ^inversely proportional to the total mass of the system [91]. For 
instance a BH-NS binary with 5M BH mass has / pea k ~ 1 — 2 kHz, depending on, 
e.g., the BH spin [76, 66]. 

This scenario is a good candidate for the creation of accretion-powered GRBs. 

• Hypermassive NS formation — NS-NS binaries with total mass < 3M (this mass 
threshold depends on the NS EOS [55, 82, 90]) will not promptly collapse into a BH, 
but will first form a so-called hypermassive NS || [95, 73, 81], supported by differential 
rotation and thermal pressure [42, 61, 90]. The hypermassive NS eventually collapses 
into a BH with a delay of 1ms- Is due to (i) losing angular momentum via GWs or 

§ Mass shedding commences much earlier, at significantly GW frequencies than tidal disruption. Note 
that Ref. [89] identified mass shedding as the cutoff frequency, which underestimates the tidal disruption 
frequency [67]. 

|| Neutron stars are called hypermassive if they exceed the mass limit of rigidly rotating NSs [92]. 
NSs that are above the mass limit of non-rotating NSs, but whose mass could be supported by rigid 
rotation, are called supramassive; e.g., [93]. For a binary with low-mass (~ 1M ) NSs, it is possible 
that the resulting post-merger NS mass can be supported even without differential rotation (i.e. it is 
not hypermassive). In this case the NS can be long lived (see, e.g., [94]). 
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magnetic processes [96] (magnetorotational instability [93] or magnetic winding [97]), 
(ii) its rotation becoming more uniform due to magnetic braking and viscosity [98, 99], 
in which case differential rotational no longer supports the NS, and/or (iii) cooling 
due to, e.g., neutrino emission, so thermal pressure provides less support against the 
gravitational pull [90, 100]. 

A rapidly rotating hypermassive NS will assume a (non-axisymmetric) ellipsoidal 
shape [81], which is energetically favorable over a spheroid [42]. Such an ellipsoidal 
hypermassive NS will emit a strong GW signal at twice its (quasiperiodic) rotational 
frequency / qpd - 2 - 4kHz [63, 55, 82, 90]. 

Hypermassive NSs may leave a massive accretion disk behind after collapsing into 
a BH [73, 82, 101]. It seems that the outcome depends on the binary mass. As 
we saw above, for high-mass NSs, the merger promptly forms a BH. At the highest 
masses for which a hypermassive NS is formed, a very short-lived NS is formed, with 
likely suppressed accretion disk formation [82]. For lower-mass binaries, the formed 
hypermassive NS is longer-lived, resulting in stronger GW emission. These lower-mass 
binaries also result in the formation of a more massive accretion disk. 
In this scenario, in which a hypermassive NS is formed, a significant amount of GW 
energy can be emitted from the hypermassive NS at around its quasiperiodic rotation 
frequency. Quasiperiodic GW emission from a hypermassive NS would be detectable 
with advanced detectors from ~ 20Mpc [102, 90], especially because it would be 
accompanied by an inspiral phase with significantly higher signal-to-noise ratio (SNR). 
This scenario is a good candidate for the creation of GRBs, with either BH-torus or 
protomagnetar central engines. 

• Formation of stable NS — If the total mass M binary of a NS-NS binary is below 
the mass limit M N s, ma x of non-rotating NSs, or if the merged NS's mass is < M N s, ma x 
due to, e.g., tidal disruption, a stable, long-lived NS can form from the merger. While 
typical observed NS-NS binary masses [103] are likely above M N s ima x ^ 2M [104], 
it is plausible that some mergers end up forming a stable NS. 

Accretion disks themselves may contribute to the GW emission of binaries [105, 
106, 107, 108]. BH-torus systems can be unstable to non-axisymmetric perturbations 
that may give rise to non-axisymmetric torus structure, resulting in GW emission at a 
few- hundred Hz. See also Section 2.2.2 for GW emission from accretion disks. 

2.1.3. Ringdown phase — When a BH is formed, or when substantial matter with 
non-axisymmetric structure (i.e., NS or fragmented accretion disk) falls into the BH, 
the event horizon of the BH is initially perturbed. Subsequently, it quickly approaches 
the non-perturbed (Kerr) solution via the emission of GWs. This process is called 
ringdown, with characteristic frequencies of / pea k ~ 6.5 — 7kHz (e.g., [55, 82]) ^f, given 

If Compare to the orbiting frequency /isco = 4.4(M/M ) _1 kHz of a particle around a Schwarzschild 
BH of mass M at the innermost stable circular orbit (ISCO). The orbital frequency greatly depends 
on the BH spin (e.g., [109] for the specific dependence). 
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Figure 3. 

collapses. 



Schematic diagram of GW emission scenarios of massive stellar core 



the typical mass of the binary (~ 2 — 3M ). 

At and above / pea k, the GW spectrum is qualitatively independent of the properties 
of the binary (although quantitatively the characteristic frequency decreases with 
increasing binary mass). It decays exponentially due to the m = 2 quasinormal mode 
oscillation of the final black hole [62, 42]. For solar-mass binaries, GW radiation from 
the BH ringdown is undetectable with advanced GW interferometers due to its high 
frequency. 

We note that BH ringdown can be suppressed or even choked in the presence of 
intense accretion onto the BH, as is likely the case for unequal-mass NS-NS binaries 



2.2. Core Collapse 

Massive stars develop an inert iron core supported by non-thermal degeneracy pressure 
in their center as the final stage of nuclear fusion. The growing iron core becomes 
gravitationally unstable upon reaching a mass around the Chandrasekhar mass (~ 
1.44 M ) due to a softening of the EOS as the degenerate fermions become relativistic 
[110], collapsing into a NS or a BH. Such a collapsing stellar system may result in 
a supernova explosion (a core-collapse supernova), and/or a GRB. The central engine 
driving GRB emission may be an accreting BH (e.g., [14]) or a rapidly spinning, strongly 
magnetized protoneutron star (a millisecond proto-magnetar] e.g., [Ill] and references 
therein) . 

The core collapse of a massive star may emit GWs through various mechanisms 
[38, 33]. Below we focus on the mechanisms producing GWs that may be sufficiently 
strong to be detected at distances relevant to GRBs, i.e. that may be detectable from 
3> 10 Mpc with 2 nd or 3 rd generation GW interferometers. A schematic diagram of these 
emission processes are shown in Fig. 3. 

2.2.1. Rotational instabilities in protoneutron stars — For massive stars with initial 
stellar masses 10 M < M < 25 M , the collapsing core is expected to form a so- 



[83]. 
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called protoneutron star [112]. The resulting rotating protoneutron star can be unstable 
to non-axisymmetric deformations, potentially giving rise to copious GW emission 
[113, 114, 115, 116, 117, 118, 119]. The onset of rotational instabilities depends on the 
rotational rate of the star, which can be conveniently parameterized by /3 = T TOt /\W\, 
i.e. the ratio of the star's rotational kinetic (T rot ) and gravitational potential (W) energy 
[120]. The resulting non-axisymmetric structure may be a bar-like m = 2 mode, giving 
rise to a characteristic GW emission. Higher m modes may also arise, albeit they have 
longer growth time [121]. While there are still various uncertainties in the evolution 
and role of rotational instabilities (e.g., non-linear mode-coupling effects can severely 
limit the deformation [117], or the role of viscosity and neutrino cooling), the emerging 
picture is that rotational instabilities are likely viable emitters of GWs and could play an 
important role in the future detection and understanding of GRB progenitors through 
their GW signature. 

The energy potentially available for GW emission is abundant. The rotational 
energy of a typical NS with 1kHz rotational frequency is ^ lO _2 M c 2 (e.g., [122]). 
Even a fraction of this energy, if radiated away in GWs, could be detectable at large 
distances (3> lOMpc) with advanced detectors. The protoneutron star may also accrete 
supernova fallback material. Such an accretion further increases the angular momentum 
and energy available that may be radiated away via GWs [123, 50, 124]. 

The amplitude of a GW signal emitted by a rotating bar scales as h ~ MR 2 f 2 /d, 
where M, i?, / and d are the mass, radius, GW frequency (i.e. twice the rotational 
frequency), and distance of the NS, respectively [50]. The energy radiated away in 
GWs, in the Newtonian quadrupole approximation, can be estimated as (see, e.g., [50]) 




where e is the ellipticity of the bar and r is the duration of the presence of the instability. 
Protoneutron stars are subject to different rotational instabilities, the two main 
categories being dynamical and secular instabilities. A star is dynamically unstable 
if it is unstable to non-axisymmetric perturbations even in the absence of dissipation 
(i.e. if a slightly non-axisymmetric shape that conserves angular momentum is 
energetically favorable). A star is secularly unstable if it is unstable to non-axisymmetric 
perturbations only if dissipative effects are relevant, i.e. if the change towards a non- 
axisymmetric shape requires the radiation or redistribution of angular momentum (see 
[125, 126] for interesting historical overviews). Below we discuss these instabilities 
further in detail. 

• Dynamical instabilities - Rapidly rotating stars will be subject to dynamical 
instabilities driven by hydrodynamical and gravitational effects [113, 127, 117]. 
Dynamical instability is the "simplest" form of NS instabilities, since its development 
is quick (on the time scale of the rotational period of the NS), and does not require 
dissipation. A uniformly rotating, classical fluid body becomes dynamically unstable 
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at a rotation rate /3 > 0.27 = /3d yn + [128]. Stability conditions are essentially 
the same for relativistic stars, for which /5d yn ^ 0.24 [117, 50]. Furthermore, 
differentially rotating stars are subject to non-axisymmetric instabilities even at 
much slower rotation with /3 < 0.09 [129, 130, 131, 132]. Such low-/3 instabilities in 
differentially rotating protoneutron stars are probably analogous to the Papaloizou- 
Pringle instability ([133]; also see Section 2.2.2): the protoneutron star core is 
surrounded by a fluid rotating at the frequency of a non-axisymmetric mode of the 
core, hence exciting this mode [127]. As numerical simulations so far have been 
too short to capture the long-term behavior of some dynamical instabilities, they 
are not conclusive in terms of the total energy emitted via GWs from dynamical 
instabilities. Nevertheless, GW emission seems to be fast relative to the cooling time 
of the protoneutron star, or compared to energy loss due to viscosity. Consequently, if 
competing mechanisms that radiate away angular momentum (e.g., magnetic fields) 
are weak, GWs can carry away a significant fraction of the protoneutron star's 
rotational energy, producing a signal that may be detectable from 

for narrow-band, circularly polarized GW signals from optimal source direction*. The 
frequency scaling of the distance only applies to / > 300 Hz. 

• Secular instabilities - Stars with lower rotation rates (i.e. for which the faster- 
developing dynamical instabilities are not prevalent) can be subject to secular (i.e. 
dissipation-driven) non-axisymmetric instabilities [135, 136, 137]. Dissipation can 
occur via gravitational radiation [136] or fluid viscosity [113, 50]. GW emission drives 
frame-dragging (so-called Chandrasekhar-Friedman-Schutz) instabilities of modes 
that are retrograde with respect to the star but prograde with respect to the observer 
[138, 139]. Among GW-driven instabilities, fundamental /-mode bar instabilities 
have the shortest growth time: 0.1s < r GW < 7 x 10 4 s for 0.27 > ft > 0.15 [135]. 
In the uniformly rotating approximation of a relativistic star, the protoneutron star 
becomes unstable to GW-driven m = 2 /-mode instabilities for f3 > 0.06 — 0.09 = /5 sec 
depending on the EOS and stellar mass [140] (compare with the Newtonian limit 
of /3 > 0.14 [138, 135]). At lower rotation rates, stars become unstable to higher- 
multipole / modes, albeit higher modes have longer growth times [121, 141]. A recent 
work of Passamonti et al. [141] indicates that the most unstable / modes developing 
in the aftermath of a supernova explosion may in fact be / = m = 3 and 4 modes. 
As angular momentum is radiated away through GWs (and potentially other 
channels), the protoneutron star's rotation frequency and therefore the emitted 
GW frequency decreases, sweeping towards the most sensitive band of LIGO-like 
detectors. Some analytical and numerical results indicate that such a GW signal 

+ Note that a uniformly rotating, constant density star cannot reach such high /?. The maximum 
rotation of a compact star is limited by mass-shedding to f3 ~ 0.1 [116]. 

* Given 10 x sensitivity improvement for advanced detectors; based on the results of [134] that used 
only one GW detector - multiple detectors could see even farther. 
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may be detectable at large distances, up to ~ 100 Mpc with advanced interferometers 
[135, 142, 143, 144]. Some other recent simulations are also promising, even though 
they only cover the first few milliseconds after core collapse and therefore cannot 
capture the long-term evolution of the bar mode [130, 131]. Nevertheless, some recent, 
realistic simulations (e.g., [141, 145, 146, 147]) predict much smaller detectable range. 
Rotating protoneutron stars are also unstable to GW-driven r-mode oscillations (the 
restoring force being the Coriolis force) at any (i.e. arbitrarily low) rotation rate 
[148, 149, 50, 150]. R modes are important only if their growth time is less than the 
damping time of viscous forces. Further, r-mode instability is expected to be saturated 
at low amplitude due to dissipative effects [65], and is suppressed in the presence of 
magnetic fields [151]. Nevertheless, under favorable conditions, GW signals from 
protoneutron star r modes may be detectable for several years after core collapse 
[152]. Given such a long duration, the emitted GW signal may be integrated for a 
measurement of ~ 1 yr that would give a detectable signal to distances of ~ 30 Mpc 
[152] (or to even farther if the protoneutron star is a strange quark star [65, 153]). 
Further calculations (see [50] and references therein) suggest that the saturation (i.e. 
maximum) amplitude of r-mode instabilities is limited by its non-linear coupling to 
other inert ial modes. 

• Magnetic distortion - Toroidal magnetic fields (B > 10 12 G) can distort a rotating 
NS into a prolate shape [154, 155, 156, 157, 158]. Such configuration is unstable to 
the growth of the angle between the NS's angular momentum and the magnetic axis. 
This angle increases until the angular momentum and magnetic axis are orthogonal. 
The resulting configuration is a rotating, non-axisymmetric body that is an efficient 
emitter of GWs [155, 159]. Fast rotating magnetars can also lose spin energy 
via magnetic dipole radiation and/or magnetized, relativistic winds. Nevertheless, 
observations indicate that such energy losses are not typical ([159] and references 
therein), suggesting that GW emission may be relevant in the early phase of newly 
born magnetars. Given that magnetic distortions indeed result in efficient GW 
emission, Dall'Osso et al. [159] find that such a GW signal can be detected out 
to the Virgo cluster with advanced detectors. 

Besides their initial rotation, protoneutron stars can accrete material from the 
infalling matter after core collapse [160]. It is possible that NSs radiate away this 
angular momentum that they gain from accretion via GWs, e.g., via dynamical and/or 
secular instabilities discussed above. Such GW signal may be detectable to > 10 Mpc 
with advanced detectors [124]. 

2.2.2. Accretion disk non-axisymmetric instabilities - Upon the core collapse of a 
massive star (> 30 M [161]), a plausible scenario is the formation of a central BH 
surrounded by an accretion disk [162, 163, 164]. Such a BH-torus system can be 
the source of copious GW emission if the disk assumes a finite quadrupole structure 
due to non-axisymmetric instabilities [38]. The emergence of such a non-axisymmetric 
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structure on the time scales comparable to the lifetime of the accretion disk requires a 
stellar progenitor with sufficiently high angular momentum. A high rotation rate is also 
required for the creation of a GRB (e.g., [165]). Below we outline some of the possible 
scenarios through which non-axisymmetric instabilities in accretion disks may result in 
strong GW emission, and address whether existing or future GW detectors can probe 
that signal. 

As an example to the sensitivity of a GW search to accretion-disk non-axisymmetric 
instabilities, we consider a 3M BH with a 0.01 M clump in its accretion disk. Such a 
system emits GWs similar to an isolated low-mass binary coalescence, and may be seen 
to a horizon distance (for definition see Section 4.1.1) [52] 

D h ~ 55 Mpc (Mc/OA M ) 5/6 (3) 

for an advanced detector with optimal source position and orientation (this distance 
is reduced by a factor 2.26 for average sky location and orientation), where A4 C = 
(mim 2 ) 3 / 5 /(mi + tt^) 1 / 5 is the so-called chirp mass (A4 C = 0.1 M for the example 
above) . 

• Disk fragmentation via gravitational instability Accretion disks with 
sufficiently large angular momenta are gravitationally unstable [166, 167, 168]. A 
gravitationally unstable disk will fragment if the disk cooling time is sufficiently 
short (< orbital period) [169]. Relevant accretion disks cool rapidly, e.g., via 
neutrinos or strong winds [34], favoring fragmentation, although the emergence of 
magnetorotational instability (MRI) in the disk can heat the disk through ohmic 
dissipation [170]. The resulting fragmented disk will emit a strong, chirp-like GW 
signal [171]. Viscosity and GW emission drive the angular momentum loss of the 
disk. Consequently, low disk viscosity favors stronger GW emission, which may be 
detectable from ~ 100 Mpc with Advanced LIGO/ Virgo [167]. The GW frequency at 
its highest SNR is probably in the 10 2 — 10 3 Hz range, depending on the strength of 
viscous forces. The duration of GW emission may be similar to that of gamma-ray 
emission. Furthermore, the tidal disruption of fragments in the accretion disk may 
be behind some of the X-ray flares observed in GRB afterglows [167]. Accretion disk 
fragmentation has been observed in approximate numerical simulations [168, 172]. It 
is not clear, however, whether such a fragmentation would occur in more realistic 
cases as well, unless the pressure support of the star is drastically reduced, e.g., via 
neutrino losses [172]. 

• Non-axisymmetric structure via Papaloizou-Pringle instability Differen- 
tially rotating accretion disks can be subject to global non-axisymmetric instabilities. 
One such instability, discovered by Papaloizou & Pringle [133], develops in accretion 
disks in which azimuthal pressure gradients (due to high internal temperature) give 
rise to differential rotation [133, 173, 174, 175], i.e. the angular velocity Q(R) of 
the disk depends on the radius R. In such disks, large-scale spiral pressure waves 
can emerge with fixed pattern speed Q p . If there is a radius R C) the so-called coro- 
tation radius, at which the pattern speed is equal to the disk rotational rate (i.e. 
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fi(i? c ) = the wave within i? c will propagate slower than the rotation of the disk, 
in fact decreasing the mechanical energy of the disk, resulting in a so-called negative 
energy wave (e.g., [175]). Such negative energy wave can interact, at the corotation 
radius, with the positive energy wave that develops at R > R c . The negative energy 
wave can increase its amplitude by losing energy, which in turn can increase the am- 
plitude of the positive energy wave, thus feeding the instability. See, e.g., [175] for an 
expressive description of the phenomenon. 

The Papaloizou-Pringle instability gives rise to an (m = 1) non-axisymmetric 
structure on a dynamical time scale (i.e. over a time period comparable to the 
rotation period). Such non-axisymmetric structure can persist for much longer than 
the dynamical time scale, resulting in strong GW emission. The Papaloizou-Pringle 
instability and the resulting non-axisymmetric structure have been observed in 3D 
relativistic simulations of BH-torus systems [176, 108]. These simulations indicate 
that BH-accretion disk systems subject to the Papaloizou-Pringle instability emit 
GWs in the 10 2 — 10 3 Hz frequency range that may be detectable up to ~ 100 Mpc 
with Advanced LIGO/Virgo [108]. Strong magnetic fields present in the accretion disk 
can enhance the instability for thick disks (and may suppress it for thin disks) [177]. 
Nevertheless, we note that the Papaloizou-Pringle instability was found in simulations 
for initially axisymmetric tori [108]. Such complete axial symmetry may not develop 
in compact binary mergers. Numerical simulations of compact binary mergers have 
not yet shown signs of the development of the Papaloizou-Pringle instability (e.g., 
[97]). 

• Suspended accretion - In order for accretion around a central, rotating BH to 
continue for the time scales of long GRBs, it has been suggested that accretion may 
be "suspended" (i.e. slowed down or temporarily halted; e.g., [178, 178, 179]). Such 
suspended accretion would be achieved through magnetic fields, which transfer some of 
the rotational energy of the BH to the disk. A fraction of this rotational energy could 
then be radiated away through GWs [180, 181]. In order to emit GWs, the accretion 
disk needs to lose its axial symmetry, e.g., through fragmentation. GW radiation from 
suspended accretion has been suggested to carry away as much as Eqw ~ 10 -2 M c 2 
in the sensitive frequency band of LIGO/ Virgo-like interferometers over a duration 
comparable to the duration of long GRBs (~ 30 s) [178, 179]. Suspended accretion, 
nevertheless, requires highly ordered accretion, which may be difficult to achieve due 
to the development of MRI heating or disk turbulence (e.g., [182]). Simulations to 
date (e.g., [Ill]) have not provided support for suspended accretion, even in the 
presence of ~ 10 15 G magnetic fields. 

2.2.3. Fragmentation of collapsing core — In very rapidly rotating stars, inf ailing 
matter may fragment even before the formation of a BH-torus system. It has been 
suggested that, in analogy with observed phenomena in star formation, the collapsing 
core may fragment and form two or more compact objects [183, 184, 171]. Such 
fragmentation was observed in relativistic numerical simulations of approximate pre- 
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supernova cores [185, 186, 187]. Such core fragmentation would give rise to strong, 
characteristic GW emission, similarly to the case of binary mergers. Nevertheless, the 
rotation rate necessary for such a fragmentation seems to be difficult to achieve with 
current stellar models [38] . 

2.3. Magnetars 

Highly magnetized neutron stars (magnetars) are thought to be the engine behind soft 
gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) [17, 188, 189, 190]. Some 
magnetars occasionally produce so-called giant flares that resemble short GRBs ([32] 
and references therein). Giant flares are much less energetic than typical short GRBs 
that are at cosmological distances, hence they are only detectable from the Milky Way or 
nearby galaxies, up to < 40Mpc with current instruments [191, 32]. The lack of excess 
short-GRB population from the direction of these galaxies indicates that magnetars may 
only be responsible for a small sub-population of short GRBs [191]. 

Giant flares are much more common (with much weaker gamma emission) than 
other types of GRBs, as indicated by the rate of observed galactic and nearby giant 
flares. Due to their high rate, they may be detectable with GW observatories even 
though they are much weaker sources of GWs than other, extragalactic GRBs. 

Giant flares (as well as SGR activity or abrupt changes in the NS spin period; [192]) 
are thought to be the result of so-called starquakes: the tectonic activity (cracking) of the 
NS crust [192], which is followed by the reconfiguration of the NS's magnetic fields [191]. 
Starquakes induce seismic vibrations in the NS, causing the observed quasi-periodic 
oscillations (QPOs) in the X-ray tails of giant flares [193, 194, 195, 196, 197, 194]. 

Alternatively, it has been suggested [198, 199] that QPOs cannot be driven by 
oscillations in the NS crust due to the quick (<~ 1 s) dissipation of the vibration via 
Alfven waves into the neutron star interior. Nevertheless, the similarities between the 
statistical properties of SGR events and earthquakes may provide further evidence for 
the crustal origin of SGR events [191]. 

NS seismic vibrations result in the emission of GWs [200, 201, 154, 202, 203], 
although this emission may be weak [204]. In the optimal scenario for GW emission, 
magnetic reconfiguration can liberate 10 48 — 10 49 erg of crustal elastic energy, which is 
the upper limit for the energy radiated away via GWs [203, 154]. This energy can be 
even higher if the NS is of strange quark matter [203] . 

Quasi-periodic oscillations (QPOs) with various frequencies from ~ 10 — 10 3 Hz and 
durations up to ~ 100 s have been observed in every giant flare X-ray afterglow so far 
[193, 194, 195, 196, 197, 194]. Some of these QPOs may be connected to the seismic 
modes of the NSs [193, 197]. Of special interest are QPOs around ^ 100 Hz, which 
fall into the most sensitive frequency band of LIGO/ Virgo. Possible GW emission 
of ~ 10 44 erg at this frequency would be detectable out to ^ lOkpc with Advanced 
LIGO/Virgo (given 10 x sensitivity improvement compared to initial detectors) [134]. 
It is possible that GWs are emitted for a significantly longer time (days to months) than 
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the observed electromagnetic QPO [205], making longer-term GW searches necessary. 

A fraction of the energy from magnetic reconfiguration may excite NS fundamental 
quadrupolar fluid modes (/ modes] e.g., [206, 207], although see [145]). If NS / modes 
are excited, they are damped by GWs on a very short time scale (~ 200 ms; [203] and 
references therein), which is shorter than most other potential damping mechanisms. 
Therefore, most of the energy in the / modes may be radiated away via GWs [203]. 
NS / modes oscillate at 1 — 2 kHz (for stiff EOS (J), not too far from the most sensitive 
band of Advanced LIGO/Virgo [201]. Given the most extreme case in which 10 49 erg 
being transferred into / modes that oscillate at ~ 1kHz, the resulting GW could be 
detected up to < 2 Mpc with Advanced LIGO/Virgo (given 10 x sensitivity improvement 
compared to initial detectors) [134]. "Stacking" the GW signals from multiple events 
can extend this distance even farther [208, 209]. If the excitation energy is comparable 
to the observed giant flare energy (~ 10 46 erg) and if not all energy is radiated away 
through GWs, only galactic sources would be detectable with advanced GW detectors. 
Furthermore, recent magnetohydrodynamic simulations of magnetized NSs indicate that 
only a small fraction of the released energy is converted to / modes, which would make 
the detection of /-mode GWs less likely [146, 147]. 

2.4- Millisecond protomagnetars 

While the commonly favored model for the central engine behind GRBs is a rapidly 
accreting BH [14, 210, 211, 212], there are observational indications [213, 214] that 
a fraction of short GRBs may originate from highly magnetized, rapidly spinning 
protoneutron stars, so called millisecond protomagnetars [68, 19, 20, 69, 70]. Further, 
outflows driven by such protomagnetars could also lead to the creation of long GRBs 
[215, 216, 111]. The emergence of a millisecond protomagnetar from various cosmic 
events is therefore a plausible step in the creation of some GRBs. 

The possibility that some GRBs are driven by rapidly spinning protomagnetars 
opens the door to alternative GRB progenitors. While traditionally considered 
progenitors, i.e. (i) the core collapse of massive stars [215, 216, 69, 111] and (ii) NS-NS 
mergers, can themselves lead to the creation of millisecond protomagnetars, other cosmic 
events, such as (hi) the accretion-induced collapse of white dwarfs [217, 218, 19, 20], or 
the merger of two white dwarfs [219, 19], are also plausible progenitors. 

Millisecond protomagnetars may emit a sizable fraction of their rotational energy 
via GWs [155, 220, 159, 221], although they may also lose angular momentum via 
magnetic processes (e.g., [93, 222]). Fast rotation can result in the emergence of 
rotational instabilities, especially if the protoneutron star is differentially rotating (see 
also Section 2.2.1). If a millisecond magnet ar indeed drives a GRB, it needs to maintain 
fast rotation for durations comparable to the duration of the GRB, potentially allowing 
sufficient time for rotational instabilities to emerge and result in GW emission. If a 
significant fraction of the rotational energy (<~ lO _2 M c 2 ; e.g., [122]) is converted into 

(J The density of NSs with stiff EOS changes relatively slowly with pressure (as opposed to soft EOS). 
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gravitational radiation, the resulting GW signal could be detected to a distance of 
~ 100 Mpc with advanced detectors (see also Section 2.2.1). 

The rotation rate of millisecond magnetars, and therefore the emitted GW strength 
and frequency, likely depends on the cosmic event they originate from. Accretion- 
induced collapse of white dwarfs that accrete mass and angular momentum from a non- 
degenerate companion may produce the most rapidly rotating protomagnetars. Due to 
continued accretion, such white dwarfs will be rapidly rotating upon collapse, resulting 
in a rapidly rotating protomagnetar, and a strong GW signal. On the other hand, 
the remnant of a binary white dwarf merger will likely collapse to a magnetar with a 
significant delay (~ 10 6 yr), allowing the remnant to lose most of its angular momentum 
via stellar winds [223, 224, 225]. Binary white dwarf mergers are therefore not likely to 
be significant sources of GWs. 

3. GRB astrophysics with gravitational waves 

This section addresses the question: What can the detection of GW signals add to our 
understanding of the physics of GRBs and their progenitors? Detectable electromagnetic 
radiation, our main source of information on GRBs, is observable only from relatively 
large distances (> 10 12 — 10 13 cm) from the central engine. On the other hand, GWs are 
created right at the central engine, and can convey information about it without being 
distorted or absorbed by matter on their way to the observer. Here we outline some of 
the questions of interest that could be addressed through observing the GW signature 
of GRBs. 

3.1. Compact binary coalescence 

• Progenitors of short GRBs - As compact binary mergers will be detectable from 
hundreds of Mpc with the Advanced LIGO- Virgo network [46] , within these distances 
it will be possible to confirm or rule out compact binaries as GRB progenitors. While 
identification based on electromagnetic signals is not always straightforward [49] , the 
detection of GWs will be likely essential to unambiguously reveal the nature of the 
progenitor [39]. By connecting each measured source with the presence or absence 
of detectable GW emission, and by exploring the connections between observed 
electromagnetic counterparts and their binary progenitors, the central engine will 
be probed in unprecedented ways. For example, observing GRBs in coincidence with 
those binaries for which numerical simulations predict the formation of accretion disks 
would indicate a connection between accretion disks and GRBs. Already, analyses 
of GW data from initial detectors targeting gamma-ray bursts have helped rule out 
binary mergers as the central engines of selected events [226, 227]. 

• Population prospects - The populations of short GRBs [228], as well as compact 
binaries [46], are highly uncertain. The observation of GWs from compact binary 
coalescences could be an effective way of determining the source population. Due 
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to the very large field-of-view of GW detectors and the very weak beaming of GW 
emission, they can be used efficiently to locate practically all sources within their well- 
determined horizon distance of hundreds of megaparsecs. As source distances can also 
be determined using GWs (see next subsection), the binary coalescence population 
could also be mapped as a function of redshift, although only with the sensitivity of 
third-generation GW detectors [229] . 

• Cosmological parameters - The detection of a binary coalescence GW signal from 
an observed short GRB could be used to determine the distance and redshift of 
the source [26, 230, 231]. In fact, the determination can be done independently of 
the cosmological distance ladder. Distances of short GRBs can be reconstructed 
to a - 10 - 30% precision for < 500 Mpc for NS-NS, and for < 1.5 Gpc for NS- 
BH binaries using the Advanced LIGO- Virgo detector network [231]. Reconstructed 
source distances can then be used to accurately reconstruct the luminosity of short 
GRBs, as well as to test Hubble relation [26], with one year of observation with 
advanced detectors potentially allowing for a ~ 2% precision [232]. Further, in 
combination with measurements of the cosmic microwave background, reconstructed 
binary merger distances could be used to constrain the dark energy equation of state 
[232]. 

• Jet angular structure - The angular structure of relativistic jets from short GRBs is 
poorly constrained [97]. The rate of binary mergers detected through GWs, together 
with the rate of short GRBs for which a binary progenitor can be confirmed with 
GWs, could be used to determine the opening angle of short GRBs. 

Further, the polarization of the inspiral GW signals in principle could be used to 
characterize the viewing angle of observed GRBs [233]. The GW polarization from 
a binary inspiral depends on the viewing angle compared to the rotational axis of 
the binary. Towards the rotational axis, the GW signal is circularly polarized, while 
the polarization becomes elliptical for off- axis observers, eventually becoming linear 
for observers in the equatorial plane. For sufficiently strong GW signals from the 
inspiral phase, reconstructing the eccentricity of the GW polarization could provide 
information on the opening angle. With a large number of measurements, the angular 
structure of the jet could be mapped as well [233]. These correlations can be compared 
to similar constraints that might follow from the detailed multiband light curve of the 
blastwave [234]. Nevertheless, due to the relatively weak change in GW polarization 
with angle, the GW polarization may only be identified for large (>> 1°) angular 
differences, and/or for very high SNR (> 100), making the utility of GW polarization 
limited. 

• Neutron Star Equation of State - Nuclear forces have a profound influence on the 
structure of NSs (e.g., [235, 122]). They determine, among other things, the relation 
between NS mass and radius, or the NS mass limit (for the rotating and non-rotating 
cases). The EOS of matter is poorly constrained at NS densities [53]. The observation 
of NSs with masses up to ~ 2 M [236] imply a stiff EOS (a NS with soft EOS would 
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not be able to support this much mass). 

The evolution of NS-NS/BH-NS mergers strongly depends on the NSs' nuclear EOS. 
Therefore, the GW signal from binary coalescences can be used to determine/constrain 
the EOS. For instance, if a NS is tidally disrupted during a binary merger, the orbital 
frequency at which tidal disruption occurs can be used to determine the radius of the 
NS. This, together with the NS mass reconstructed from the inspiral gravitational 
waveform [25], can be used to constrain the nuclear EOS [89, 237, 76, 66, 67, 238]. 
The tidal deformation of NSs in a binary system can affect the gravitational waveform 
even prior to the merger phase. The GW energy spectrum during the last few orbits 
of a NS-NS binary prior to merger can be used to determine the compactness ratio 
(mass/radius) of the NSs [239, 238]. For these last few orbits of the binary, the GW 
frequency is at a high but still reasonably sensitive frequency for LIGO-type detectors, 
making the analysis of these last few orbits feasible. 

The tidal deformation of a NS in the inspiral phase can be described by one parameter 
(the so-called Love number [53]), which is effectively the ratio of the star's induced 
quadrupole moment to the quadrupole moment of the perturbing tidal gravitational 
field of the binary companion. Flanagan & Hinderer [53] showed that the nuclear 
EOS of NSs can be constrained even through the early inspiral phase due to the effect 
of tidal deformation on the waveform. Read et al. [54] simulated the inspiral phase of 
NS-NS binaries, showing that the NS EOS can be constrained (the NS radius can be 
determined to within <~ 1 km precision, which, together with the NS mass, would rule 
out a part of the EOS parameter space) with advanced GW detectors for a source at 
a distance of 100 Mpc. BH-NS mergers can similarly be used to gain information on 
the NS EOS [240, 67]. 

For some NS-NS mergers, collapse to a BH is delayed and a hypermassive NS is 
formed (see Section 2.1.2). During the merger, non-axisymmetric oscillation modes 
of the hypermassive NS are excited, resulting in GW emission [241, 242, 243, 244, 
95, 245, 246, 247]. The waveform of the resulting gravitational radiation depends on 
the nuclear EOS. Consequently, the observed waveform can be used to significantly 
constrain the EOS [243, 244, 95, 245, 246, 247]. The GW signature of hypermassive- 
NS oscillations could be detected with SNR= 2 out to 20 — 45 Mpc (depending on the 
EOS), which may be sufficient to constrain the EOS (at such distances, the source 
can be confirmed with high accuracy via GWs from the inspiral phase) [246] . 
Further, it is possible that the inner part of a NS becomes strange quark matter at 
high densities. A star made of strange quark matter can be self-bound, a marked 
difference from hadronic NSs that are gravitationally bound [35]. The GW emission 
of a strange quark star is substantially different from that of hadronic NSs [237], 
hence GWs could be used to determine whether NS interiors may contain strange 
quark matter [248, 249, 250, 100]. 

• Magnetic fields in protoneutron stars - Magnetic fields are suspected to play a 
significant role in the origin of short GRBs [251, 42, 97, 85], and are known to vary 
between individual pulsars (e.g., from their spindown rate [252]). In the inspiral phase, 
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internal magnetic fields in NSs (for realistic field strengths) have negligible effect 
on the dynamics of the binary coalescence [253, 251, 81, 254, 78]. While magnetic 
tension in the NS can reduce tidal deformation, this effect will be too weak to be 
detectable with planned GW detectors for realistic pre-merger magnetic strengths 
(< io 14 G) [254]. For BH-NS binaries, magnetic fields seem not to alter the dynamics 
and therefore the gravitational waveform [78]. 

Magnetic fields can significantly affect the post-merger behavior of NS-NS binaries, if 
a hypermassive NS is formed after the merger. Upon merger, magnetic fields will be 
amplified via Kelvin-Helmholtz instabilities in the shear layer between the merging 
NSs [255, 81], or later via differential rotation [93, 251]. Even if magnetic fields 
are small prior to merger, the amplified fields can be strong enough to affect the 
evolution of the merger remnant and therefore GW emission. Magnetic fields in 
the hypermassive NS compete with GWs in dissipating angular momentum from 
the hypermassive NS. The decrease of angular momentum eventually results in 
the "delayed collapse" (~ 10ms— Is, depending on the mass and the EOS) of the 
hypermassive NS into a BH (e.g., [251, 246]). 

By probing the magnetic properties of a large number of NSs in binary mergers, one 
could determine, e.g., the highest achievable magnetic field frozen in a stationary NS 
([35]; e.g., via the time difference between binary merger and GRB). 

• Accretion disks The large-amplitude, quasiperiodic GW emission from 
instabilities in the accretion disk, probably together with the detection of the earlier 
inspiral phase from the binary, can be utilized to reconstruct properties of the 
accretion disk, even if no electromagnetic signal is observed from the source. Beyond 
reconstructing the mass and lifetime of the accretion disk, the development of the 
instabilities is sensitive to the properties of the accretion disk (see Section 2.2.2), i.e. 
the development of the instabilities by itself can already be informative. 

3.2. Core collapse 

• Progenitors of long GRBs - It is believed that the origin of at least some long 
GRBs are the core collapses of massive stars [256]. This could be confirmed if the 
GWs in coincidence with a long GRB progenitor were detected. Nevertheless, if no 
GW emission is detected in coincidence with long GRBs, this does not rule out the 
core-collapse model, although may constrain core-collapse dynamics. If GWs are 
detected in coincidence with a long GRB, the internal evolution of a core collapse, 
which may largely depend on the progenitor, could also be examined. For example, 
the formation of a protomagnetar or a BH-torus system, or the fragmentation of the 
massive stellar core may be differentiated via GWs. 

• Neutron star equation of state - The mass, radius and rotation rate of NSs have 
a profound effect on their potential GW emission through rotational instabilities. 
Detecting the GW signature of these instabilities can be used to constrain the mass- 
radius relation, rotation, and therefore the EOS, of NSs (e.g., [257]). 
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Further, rotational instabilities in NSs result in qualitatively different GW emission 
for conventional and strange quark NSs [65, 153]. Whether a NS is composed of 
strange quark matter could be inferred from the detected GW signature of the NS's 
rotational instability. 

• Neutron star internal physics - The evolution of the gravitational waveform from 
protoneutron star rotational instabilities strongly depends on the nuclear EOS at NS 
densities, as well as the physical parameters of the protoneutron star. Differential 
rotation inside a protoneutron star [129, 130, 131, 132], temperature (e.g., [152]), 
viscosity and neutrino cooling [152] may all leave their mark on the evolution of 
rotational instabilities and the resulting GW emission. Further, if a hypermassive 
(or supramassive; see Section 2.1.2 for their definition) NS collapses into a BH due 
to either losing angular momentum or accretion, the rotational frequency prior to 
collapse provides information on the NS EOS [237]. This rotational frequency may 
be inferred from the GW signal from the rotating NS [257, 124]. 

• Accretion physics - Because GW emission from massive stellar collapses can be 
produced by the accretion onto the central object following the collapse, the GW 
signal therefore can provide insight into accretion physics: (i) the mechanisms that 
transport angular momentum; (ii) instabilities; and (iii) interaction with the black 
hole. 

For BH-torus systems the strength of viscosity (and possibly the strength of other 
processes through which the torus loses angular momentum) can be determined, as 
these are competing effects for angular momentum loss: The total loss of angular 
momentum can be inferred from the rate of change of the GW frequency, while the 
loss through GWs can be inferred from the GW signal amplitude. 
Accretion of matter by the protoneutron star results in extended GW emission due to 
the intake of angular momentum [258, 259]. As the time scale and nature of accretion 
is likely different for accretions from supernova fallback material or material from 
a companion star, the extended GW signal from rotational instabilities will carry 
important information on the accretion mechanism. 

• Magnetic fields in neutron stars - Sufficiently strong magnetic fields inside a 
protoneutron star can be competing with GWs in radiating angular momentum away. 
As NS spindown may be reflected in the GW frequency, measuring the GW amplitude 
and NS spindown may provide information on the strength and nature of magnetic 
fields present in the NS. 

3.3. Magnetars 

The detection of GWs in coincidence with a giant flare from a magnetar could provide 
information on the processes that lead to flaring in magnetars (i.e. tectonic activity [192] 
or global reconfiguration [154] vs origin from the magnetosphere [198]). Non-detection, 
nevertheless, does not necessarily rule out these models. Upon detection, the frequency 
of GWs from NSs excited by, e.g., starquakes, could be used to infer the NS mass, radius 
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and EOS [260, 261, 262], or potentially even the strength of magnetic fields inside the 
NS [205]. 

3.4- Millisecond protomagnetars 

The detection of millisecond magnetars' GW signature in coincidence with a GRB would 
confirm their presence and role as the central engine behind some GRBs. Further, the 
GW signature of millisecond magnetars may indicate their origin (e.g., typical rotation 
rates may differ for different mechanisms; see Section 2.4). Advanced GW detectors will 
also be able to determine, out to ~ 450 Mpc, whether the progenitor of a short GRB was 
a binary merger of a NS with another NS or a BH [46] . The confirmation of a compact 
binary merger would rule out the other scenarios in which millisecond magnetars can 
form, e.g., white dwarf accretion-induced collapses, or white dwarf binary mergers. 

4. Observational strategies and prospects 

The success of detection of GWs from GRB progenitors depends on the strength of the 
emitted GWs, the sensitivity of GW observatories, as well as the observation strategies 
used to separate GW signals from the background. The search strategy, e.g., the use of 
multiple messengers, can also add to the information one can extract from the source. 

Various transient GW search strategies exist, aiming to find the GW signature of 
GRBs. These strategies differ mainly in their prior assumptions on the gravitational 
waveform. In the most general case, one can look for so-called GW bursts with 
minimal assumptions on the waveform, namely defining a maximum signal duration 
and bandwidth [263, 264, 265]. These generic burst searches can be modified to include 
additional information about the source. For instance if one considers an accretion-type 
emission, circular (on-axis) or linear (off-axis) polarizations can be required from GW 
signals, increasing detection sensitivity for signals satisfying these conditions. Other 
model-dependent assumptions include, e.g., GW signals rapidly evolving characteristic 
frequency [266], such as in the case of compact binary mergers. 

For compact binary coalescences, the gravitational waveform for a given mass 
configuration can be accurately calculated for the initial inspiral and final ringdown 
phases, and to some extent the merger phase [267]. For such systems with properties 
in a highly limited parameter space, it is beneficial to use a signal waveform template 
bank and perform matched filtering based analysis (e.g., [268]). 

Below we review the GW search strategies targeting the progenitors of GRBs. 
We indicate the prospects of such searches through presenting projected sensitivity 
estimates based on observational results with initial detectors. Advanced detectors will 
be ~ 10 x more sensitive than initial detectors, the exclusion distances therefore will be 
~ 10 x greater, corresponding to a ^1000-fold improvement in achievable source rate 
constraints. Further, advanced detectors will be sensitive in a wider frequency band. 
Additionally, the advanced detector network will eventually be larger, with the inclusion 
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of KAGRA and probably LIGO India. These improvements will further increase the 
achievable sensitivity of the advanced detector network [9], and will further improve 
signal parameter reconstruction (e.g., source direction; [269]) [9]. 

4.1. Gravitational wave search strategies 

Search strategies for GW transients can be divided into four main categories based on 
their use of other messengers. 

1 Untriggered searches, in which no information is used from other messengers. The 
analyses use GW data alone to find astrophysical sources (e.g., [270, 268]). 

2 Externally- triggered searches, in which one specifically looks for GW signals 
from sources confirmed via other messengers. For instance one can search for GW 
signals from the progenitors of detected GRBs, using their location, time, and other 
parameters (e.g., [271, 272, 273, 274, 275, 209]). 

3 Electromagnetic follow-up searches, in which GW signal candidates are used to 
trigger electromagnetic follow-up searches with other telescopes (e.g., [276, 277, 278, 
279]). 

4 Multimessenger searches, in which one uses (typically sub-threshold) signal 
candidates from multiple observatories of different messengers in a joint search (e.g., 
[280, 281, 282, 283]). 

In the following we briefly review some of the GW search strategies and past GW 
searches for GRBs, organized along the above categories. 

4.TT Untriggered searches — Untriggered transient searches provide a generic way 
to identify plausible GW sources without relying on the detection of other messengers. 
While GRBs are highly beamed (e.g., [284, 285, 286, 287, 288, 289, 290]), the beaming 
of GW emission from GRB progenitors is weak (e.g., [233]). For realistic GRB beaming 
angles, there will likely be more observations of GWs from off- axis GRB progenitors 
than GW observations triggered by GRBs [291]. 

To help quantify the source population detectable with untriggered (and other) 
searches, it is useful to characterize the sensitivity of a GW observatory to a GW source 
with the so-called horizon distance D^. The horizon distance is defined as the distance 
at which a GW source from optimal orientation (such that emission towards the Earth is 
the strongest) and optimal location (i.e. the most sensitive direction of the GW detector) 
is observed with a given SNR (typically chosen to be p = 8 for untriggered searches; see 
also Section 4.1.2) by a single GW detector. Compared to this optimal orientation and 
location, the orientation and direction-averaged distance at which the source produces 
the p = 8 is ~ Dh/2.26 (both direction and orientation- averaging contribute a factor 
- 1.5). 
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For a NS-NS binary with 1.35 M mass for each NS, the projected horizon distance 
for untriggered searches with advanced detectors (upon non-detection) is [52] 

D untrig. _ ^ Q (^/^ M ) 5/6 , (4) 

where Ai c = (mim 2 )^ 5 /(% + rn 2 ) 1 ^ is the so-called chirp mass (i.e. the horizon 
distance only depends on this specific combination of the masses of the two compact 
objects). For higher-mass binaries, the horizon distance is somewhat lower than the 
prediction of Eq. 4 (see [291]). 

4-1.2. Externally triggered searches — External triggers in GW searches provide 
additional information and increased search sensitivity for GW analyses [292, 293]. For 
example an external GRB trigger reduces the temporal and spatial extent in which one 
has to search for a GW signal [294] . A detected external trigger can also be used to set 
constraints on the GW signal. For example a short GRB implies that the GW signal is 
probably a binary merger. 

Chen & Holz [291] estimates that for GRB beaming factor < 7.5 (corresponding 
to opening angle 9 = 30°), the number of triggered and untriggered detections of GWs 
from GRB progenitors would be comparable. Further, even for greater beaming factors 
(i.e. smaller opening angles), externally triggered searches will enable the detection of 
some additional, particularly interesting, sources. The extra information in external 
triggers can be taken into account in the horizon distance by effectively lowering the 
SNR threshold (i.e. by keeping the false-alarm rate constant). The results of Chen & 
Holz imply that the sensitivity of externally triggered searches is greater than that of 
untriggered searches by a factor ^1.3 (for inspiral searches; this factor is similar for 
other types of searches as well), i.e. 

jjext.trig. _ l^ntrig^ ^ 

The same sensitivity increase was found by Dietz et al. [295]. While this increased 
horizon distance means that many more sources fall within observable reach, the actual 
number of detected sources will be decreased by the beaming of the electromagnetic (or 
other) signal, as well as the sky coverage of the available electromagnetic telescopes. On 
the other hand, all triggered sources will be detected face-on, while untriggered sources 
can take any orientation, which adds to the sensitivity of triggered searches compared 
to untriggered searches (for binary mergers, the GW amplitude is ~ 1.5 times greater 
face-on than averaged over all directions). Chen & Holz [291] find that the externally- 
triggered observation rate is a sizable addition, although less than, the untriggered GW 
observation rate, while Kelley et al. [296] find that GRB-triggered searches present no 
appreciable addition to the total rate [296]. While both of these studies approximate 
GW background as purely Gaussian, real GW background features a non-Gaussian tail 
that may increase the importance of externally triggered and follow-up searches. 

Previous GRB-triggered GW searches include searches for both unmodeled GW 
bursts and compact binary mergers. For example the short GRBs 070201 and 051103 
had electromagnetic positions overlapping nearby galaxies (Andromeda and M81, 
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respectively). No binary merger counterpart was found in either case, ruling out a 
binary progenitor in those galaxies, favoring an SGR model (a binary progenitor from 
a more distant galaxy is not ruled out, albeit it is unlikely) [297, 298]. Nevertheless, it 
is worth noting that, given their expected rate of occurrence, a binary merger within 
the distance of Andromeda or M81 (< 4Mpc) would be extremely unlikely. For GRB- 
triggered GW searches with initial LIGO/ Virgo, see also, e.g., [299]. 

Several GW searches aimed to identify GWs in coincidence with SGR flares 
[202, 134, 208, 209]. Most recently Abadie et al [209] used 1279 flares from 
six magnetars as triggers, aiming to identify GW signals from neutron-star /-mode 
ringdowns or other GW-producing mechanisms. The search set constraints on the energy 
of GW emission from these flares comparable to some giant flares' electromagnetic 
energies, and an order of magnitude below previously existing limits. For a nearby 
magnetar (SGR 0501+4516) located at ~ 1 — 2kpc that emitted a large number of 
flares during the observation period, the obtained upper limit on /-mode GW emission 
was rsj 10 47 erg. The upper limit for GW emission at ~ 10 2 Hz, i.e. within most sensitive 
frequency band of LIGO/Virgo, was ~ 3x 10 44 erg, a promising limit compared to some 
theoretical upper limits on GW emission from SGRs [203]. With Advanced LIGO/Virgo, 
one can expect, beyond more stringent limits on /-mode emission, constraints on the 
maximum emissible total energy from other NS oscillation modes (we note that these 
upper limit are constraining only if a substantial fraction of total energy goes into GWs, 
which is debated; e.g., [145]). 

4.1.3. Electromagnetic follow-up — With the rise of a global GW detector network, 
it became possible to reconstruct the direction of a GW signal candidate, albeit with a 
relatively large uncertainty. Such reconstruction enables the use of GW signal candidates 
in triggering electromagnetic follow-up observations [276, 277, 278]. There are various 
scientific advantages of such a follow-up search: 

• As electromagnetic observations cannot continuously cover the whole sky with high 
sensitivity, GWs can be used to guide these telescopes and point towards the right 
direction at the right time. 

• Often the most interesting and strongest emission from an electromagnetic transient 
occurs in a short period after its onset. Since GW emission is, in many cases, 
expected to precede the onset of electromagnetic emission, telescopes have the chance 
to commence follow-up observation in a very early emission stage, or even catch the 
onset of the electromagnetic event. 

• GW signal candidates will often have too low significance to be unambiguously 
identified as extraterrestrial signals. The observation of an electromagnetic follow- 
up event can enhance the significance of the joint observation, making detection more 
probable. 

• Similarly to other multimessenger searches, information from the different messengers 
can enhance the information (and therefore science) that one can extract from the 
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source. 



Electromagnetic follow-up searches of GW event candidates is one of the most promising 
and one of the fastest-evolving areas of GW astrophysics (see [300, 301, 302] for recent 
observational development with the LIGO- Virgo GW detectors). See Section 4.2 below 
for a more detailed description of the potentially detectable electromagnetic counterparts 
and search strategy. 

The sensitivity of electromagnetic follow-up searches will differ from the sensitivity 
of externally triggered searches due to (i) the limited reach of electromagnetic transient 
observations and (ii) the finite sky area that is scanned by observatories following up 
the GW event candidate. Taking (i) into account in defining the horizon distance, one 
gets 

D follow-u P _ ^ ^ext.trig. jj^ (g) 

where D EM is the reach of the electromagnetic follow-up observatory. Note that ]j ext - trig - 
refers to on-axis sources, while follow-up observations are the most interesting for off-axis 
sources (due to the greater source rate). For off-axis directions, GW emission is typically 
weaker by a factor ~ 1.5. The rate of follow-up observations further depends on sky 
coverage of the follow-up observatories, as well as the beaming of the electromagnetic 
signal. For a survey of the electromagnetic detectability of GW event candidates, see 
[9]. 



1^.1.^. Multimessenger analyses — Multimessenger analyses generalize the idea behind 
externally triggered searches by combining information from sub-threshold signal 
candidates of various messengers into one, more powerful analysis. While in many 
cases each type of sub-threshold messenger would individually have too low significance 
to be identified as an astrophysical signal of interest, the combination of individual 
significances can greatly enhance the analysis' potential of identifying astrophysical 
events (see Section 4.3.2 for a multimessenger search). 

The Astrophysical Multimessenger Observatory Network (AMON) initiative [303] 
plans to combine sub-threshold triggers of multiple messengers in a low-latency coherent 
multi-messenger analysis. AMON is planned to combine multiple messengers to 
identify candidates that would be sub-threshold events for individual observatories. 
The extracted information would be used, besides claiming detection itself, to initiate 
electromagnetic follow-up observations, which can add to the number of coincidentally 
observable messengers. 

The sensitivity of multimessenger searches depends on the false alarm statistics 
of the different detectors included in the search, as well as the types of messengers. 
Transient messengers with amplitude signatures, such as GWs or high-statistics photon 
signals, will have a relatively well defined horizon distance, beyond which the detection 
of such a messenger is highly unlikely. Other, discrete messengers, such as high-energy 
neutrinos, will have no such horizon distance, since there is a finite probability of 
detecting at least one neutrino even from a very distant source with expected signal 
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strength <C 1 neutrino. These distant sources can significantly contribute to the total 
number of detectable multimessenger sources, especially in the case in which their 
presence can be confirmed with other messengers (see, e.g., [283]). For such discrete 
messengers, the GW-multimessenger horizon distance is comparable to that of externally 
triggered searches: 

jjmultimessenger ^ j-^ext.trig. ^ij^ 

The rate of detected multimessenger sources will be limited by beaming, and the 
probabilistic nature of observing at least 1 discrete messenger (e.g., high-energy neutrino 
[282]). 



4-2. Electromagnetic counterparts 

This section outlines the observable electromagnetic emission of GRB progenitors, 
focusing on counterparts other than gamma rays, and their utility in GW searches. 

Due to the highly beamed emission of most GRBsff , the majority of GRBs are off- 
axis, i.e. their prompt emission cannot be observed from the earth. Nevertheless, weaker, 
off-axis electromagnetic emission from GRB progenitors could be detected during GW 
follow-up observations, if the electromagnetic emission is present beyond a few minutes 
after the prompt emission (i.e. the delay between detecting a GW signal candidate and 
pointing a telescope towards the expected source location). 

Most compact object mergers observed by LIGO /Virgo will not be accompanied 
by observable GRBs due to the small gamma-ray beaming angle. For this reason, 
the last few years have seen extensive investigation into rapid follow-up search 
strategies for more isotropic, but more difficult to identify, electromagnetic counterparts 
[277, 278, 304, 300, 301, 305]. Follow-up searches aim to minimize the time between 
the identification of a GW event candidate and the start time of the observation with 
electromagnetic telescopes in the required directions. The time delay of sending a request 
for follow-up search to electromagnetic telescopes after a GW event was ^tens of minutes 
for initial searches (mostly due to the manual confirmation of event selection) . This delay 
will likely be reduced to ^minutes for searches with advanced detectors [306] . 

The uncertainty in direction reconstruction of GW event candidates with the 
planned advanced detector network will likely be in the ^tens of square degrees range 
(e.g., [269, 9]). With such an uncertainty, telescopes with large fields of view will need 
to be utilized, focusing on electromagnetic emission that is ongoing for hours-days so 
the telescopes can scan through the interesting sky area. Existing galaxy catalogs can 
also be used to largely decrease the sky area that needs to be mapped (see Section 4.4). 
We note here that, besides GW observations, the observation of the prompt gamma- 
ray emission itself has a position resolution of a few degrees, i.e. comparable to the 

ff Beaming of short GRBs in general is not well constrained [13]. Nevertheless, beaming has been 
confirmed for a few short GRBs [287, 290], indicating a highly beamed jet with half opening angles 
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resolution of GW observations. Afterglow observations are therefore also important for 
the direction reconstruction of observed GRBs. 

Beyond direction reconstruction, GRB afterglows also promise to carry important 
information connected to GW radiation. Corsi & Meszaros [137] pointed out that 
the observed X-ray plateau in the afterglow of some GRBs may be connected to the 
formation of a highly magnetized millisecond pulsar that can lose angular momentum 
through GWs during a 10 3 — 10 4 s period (see also [307]). 

Several promising emission models have been identified which predict longer-lived, 
weakly (or not at all) beamed emission that can be detected on a distance scale 
comparable to the sensitivity range of advanced GW detectors. In the following list 
of these models, we focus on optical counterparts of compact binary coalescences. For 
a survey of the detectability of some of these emission models, see, e.g., [9]. 

• Electromagnetic remnants from compact binary mergers — NS-NS and NS- 
BH mergers can eject energetic (sub-relativistic or relativistic) outflows. These 
energetic outflows can interact with matter in the interstellar medium, creating a 
long-lasting radio signal [308, 309]. Radio remnants from dynamically ejected sub- 
relativistic material can appear months to years after the binary merger, and can last 
for several years [309] . A deep radio survey can identify these radio remnants out to 
~ 300 Mpc [308, 309]. Mildly /ultra-relativistic outflows may create even brighter 
radio emission, on the time scale of weeks [309]. The detectability of the radio 
transient strongly depends on circum-merger density, and therefore the signal may 
be lower for binaries that left their galaxy prior to merger. 

An interesting case is shocks generated within the NSs upon the merger of a NS- 
NS binary [90, 310], which can plausibly drive ultra-relativistic outflows. These shock 
waves, besides heating the NSs, result in shock breakout from the NSs' surface, driving 
a nearly omnidirectional, ultra-relativistic outflow [311]. The outflow decelerates in 
the ambient medium, resulting in a bright X-ray flare seconds after the merger that 
could be detectable with, e.g., Swift XRT [312]. The flare gradually changes to optical 
and radio as the outflow decelerates, potentially detectable with, e.g., Pan-STARRS 
[313] and EVLA [314]. 

• Macronovae/kilonovae — Matter ejected from compact binary encounters may 
produce electromagnetic optical or near-infrared transients that are suitable for follow- 
up searches of GW signal candidates [315, 316, 317, 318, 319, 320, 321, 47, 101, 309]. 
Matter ejected by binary mergers (tidal tails or accretion disk outflows; e.g., 
[316]), initially at nuclear densities, expands and undergoes r-process nucleosynthesis, 
producing heavier, radioactive elements. The decay of such elements produces an 
isotropic emission that lasts for days, named "macronova" [316] or "kilonova" [318]. 
This is an appealing model, because it suggests that virtually all LIGO detected NS- 
NS or BH-NS mergers would have detectable optical counterparts. However, finding 
these counterparts in the large GW directional uncertainty region would likely require 
a combination of powerful survey telescopes and follow-up spectroscopy, as used by 
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Palomar Transient Factory, and Pan-Starrs, or the future LSST [322, 296, 309, 47]. 

• X-ray /optical afterglow — The relativistic outflow from a GRB central 
engine interacts with matter in the circumburst medium, emitting a long-lasting 
electromagnetic afterglow in a wide range of frequencies (e.g., [234, 47]). At X- 
ray and optical frequencies (lasting from minutes to a day) , afterglow emission is the 
strongest shortly after the prompt GRB in the direction of the relativistic outflow (i.e. 
on- axis), and rapidly decays afterwards [234]. On-axis X-ray emission is detectable 
(e.g., with the Swift XRT telescope) over a period of ^ 1 day after the prompt burst 
for source distances out to the horizon distance of GW telescopes [305] . Off- axis X-ray 
and optical emission, far from the axis of rotation, peaks at a significantly lower flux, 
making the detection of off- axis X-ray less useful in multimessenger searches [234]. 
For radio frequencies around ~ 10 GHz, on-axis emission for long GRBs peaks at 3-6 
days in the source rest frame after the prompt burst [323], while off- axis emission 
peaks on a similar time scale. 

• Magnetic interaction between BH & NS — For NS-BH binaries, once the BH 
enters the magnetosphere of the NS, it can interact with the NS magnetic field [84]. 
For strong enough NS magnetic fields (~ 10 12 G), such an interaction results in the 
copious emission of electromagnetic radiation, via a mechanism analogous to the 
Blandford-Znajek mechanism [324]. The electromagnetic radiation may resemble an 
extremely short (C?(ms)) GRB , and could be detected from distances of > 60Mpc 
using the Swift or Fermi satellites [84]. Although this mechanism produces a short 
electromagnetic emission, it is expected to be less beamed than the "standard" 
gamma-ray emission. 

• Magnetic interaction between NS & NS — The magnetic interaction of NSs in 
NS-NS binaries can drive an electromagnetic (X-rays or gamma-rays) signal within 
the last few seconds prior to merger through electric dissipation [85]. As the non- 
magnetic NS moves through the magnetosphere of the other, highly magnetized NS 
(~ 10 12 — 10 14 G), the induced electromagnetic force across the non-magnetic NS sets 
up a circuit connecting the two stars. Some of this electric potential dissipates either in 
the surface layer of the magnetic NS or in the space between the two stars, depending 
on the resistivity of the space between the two stars. A part of the dissipated energy is 
carried away through electromagnetic emission, which can reach up to <~ 10 44 ergs _1 
in the last < Is prior to merger [325, 326]. 

• Neutron star flares induced by tidal crust cracking — NSs in inspiraling 
binary systems are subject to tidal deformation. When this deformation exceeds 
a critical level, the crust of the NS cracks, possibly resulting in a NS flare similar to 
observed flares [327, 328, 329, 330] (see Section 2.3 on magnetar giant flares). Such 
flare would represent an isotropic electromagnetic counterpart for binary mergers. 
Observed short-GRB precursors may be created by crust cracking [329]. 
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4-3. Neutrino counterpart 

Many GW sources, and in particular GRBs, are expected to be copious emitters of 
neutrinos [331, 332, 333, 334, 335, 336]. Astrophysical neutrino emission is expected 
within two main sub-groups based on two distinct emission mechanisms, with two 
distinct energy ranges. MeV energy neutrinos (with energies t v < 100 MeV) are 
produced in the extremely hot, dense central regions of core-collapse supernovae and 
probably GRBs. High-energy neutrinos (with energies e u > 100 GeV) are expected to 
be emitted by shock accelerated particles in relativistic outflows driven by the central 
engine of the GRB (e.g., [333]). So far only MeV energy astrophysical neutrinos have 
been confirmed, and in one instance, for supernova 1987A [337, 338]. 

One of the advantages of joint GW-neutrino searches [282, 283] is that GW and 
neutrino detectors (MeV and high energy) continuously observe the whole sky (one high- 
energy neutrino detector observing half the sky) , recording signal candidates without the 
need to "point" the detector in a particular direction. Such full sky coverage is of special 
importance for multimessenger searches as the sky coverage for each messenger needs 
to overlap for a joint search. Further, while electromagnetic follow-up searches require 
low-latency response from electromagnetic telescopes, GW-neutrino observations can be 
performed with high latency without loss of information. Nevertheless, fast analysis of 
joint GW-neutrino observations can also enable the electromagnetic follow-up of joint 
event candidates. 

4.3.1. MeV neutrinos — In core-collapse supernovae, most of the released gravitational 
binding energy is emitted in a burst of ^MeV neutrinos (e.g., [339]). Both MeV neutrino 
[340, 341] and GW emissions [33, 342] are expected to commence near core bounce within 
a few milliseconds (< 10 ms). This temporal correlation is orders of magnitude tighter 
than correlation with electromagnetic signals [281, 336], and can greatly enhance the 
sensitivity of a search for joint GW-MeV neutrino sources. MeV neutrino emission peaks 
within a fraction of a second after bounce, while emission continues for up to tens of 
seconds as the protoneutron star cools and contracts [341, 339]. 

The joint detection of GWs and MeV neutrinos could provide constraints on the 
core-collapse supernova mechanism as well as information on the properties of matter 
at high energies and densities [336]. For example the neutrino spectrum from core- 
collapse supernovae depend both on the nuclear EOS [343] as well as the spin of the 
core [344] . One can break this degeneracy using the additional information available in 
the GW channel, inferring information on both the EOS and the spin of the core. A 
potential secondary collapse due to hadron-quark phase transition in the protoneutron 
star during its post-bounce evolution [345] would also result in characteristic neutrino 
and GW emissions [345, 346]. 

Long and short GRB central engines are also probably strong emitters of MeV 
neutrinos [331, 332, 163, 164, 164]. As the progenitors of long GRBs are likely the 
progenitors of core-collapse supernovae as well [16, 256, 347], emission from long GRBs 
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may commence similarly to the case of core-collapse supernovae. In the core-collapse 
scenario, neutrino emission from the core ceases upon black hole formation, while a 
newly formed accretion disk will become a significant source, as a large fraction of the 
accretion energy can be carried away by neutrinos [331]. 

Compact binary mergers - the likely progenitor of short GRBs - are expected to 
emit MeV neutrinos [331, 332] from the accretion disk around the central black hole. In 
the case of NS-NS binary mergers, a hypermassive NS may form prior to BH formation, 
and further produce neutrino [332] as well as GWs [81] emission. MeV neutrino emission 
from compact binary mergers is expected to commence at the beginning with the merger 
phase, in coincidence with the GW burst due to the merger. This GW burst signal is 
preceded by the much longer inspiral phase, and followed by the ringdown phase. Such 
temporal coincidence can be utilized in a joint search. Nevertheless, for typical distances 
of binary mergers, MeV neutrinos would be difficult to detect [348]. 

Several large-scale MeV neutrino detectors are currently in operation. These 
include Super-Kamiokande (Japan) [349], KamLAND (Japan) [350], LVD (Italy) [351], 
Borexino (Italy) [352] and Baksan (Russia) [353]. Further, the IceCube high-energy 
neutrino detector [354] is also capable of detecting bursts of MeV neutrinos [341], albeit 
without the ability to reconstruct the source direction. Recognizing the importance 
of multimessenger observations, Super-Kamiokande, LVD, Borexino and IceCube are 
members of the Supernova Early Warning System (SNEWS) [355]. These observatories 
send real-time triggers of detected supernova candidate events, which are distributed to 
the astronomer community to allow for low-latency follow-up electromagnetic (or other) 
searches. 

The most sensitive current MeV neutrino detectors (Super-Kamiokande and 
IceCube) are able to detect the MeV neutrino signal of a supernova from up to 
^100 kpc [356]. The expected event number within this distance, however, is relatively 
small outside of the Milky Way. Planned megaton detectors, such as LBNE [357] or 
Hyper-Kamiokande [358], could detect supernovae from up to <10 Mpc [359] with 
a supernova rate of ^ 1 per year. Multimessenger searches of MeV neutrinos and 
GWs could further increase these detectable supernova rates [360], and would provide 
increased confidence in a detected signal. This can be especially important when no 
electromagnetic counterpart is detected. 

4.3.2. High-energy neutrinos — Non-thermal, high-energy (3>GeV) neutrinos are 
thought to be created within relativist ic outflows driven by the central engine [333, 
361, 334, 362, 363, 364, 365, 366, 367]. Emission from early GRB afterglows is also 
plausible (e.g., [368]). The emission mechanism is likely similar for both long and short 
GRBs, the former expected to be the stronger emitter. Core-collapse supernovae with 
rapidly rotating cores [366, 369], magnetars [370, 371] and millisecond protomagnetars 
[372] are also thought to emit high-energy neutrinos. While recent limits obtained with 
the IceCube detector constrain some emission models [373], the standard fireball picture 
for neutrino emission remains viable [374, 375]. 
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As their production mechanism is connected to that of gamma photons, high-energy 
neutrinos are probably beamed similarly to that of gamma ray emission. This decreases 
the rate of detectable joint sources similarly to the case of electromagnetic emission. 

For a joint detection, an especially interesting class of sources are GRBs with little 
or no detectable electromagnetic counterpart, which are therefore only detectable via 
GWs and neutrinos. Such sources include so-called choked GRBs. In the core-collapse 
supernova scenario, gamma-ray emission is observable only if the relativistic outflow 
from the central engine, that is responsible for the production of gamma-rays, breaks 
out of the star [376, 283]. The outflow, in order to advance, needs to be driven by an 
(active) central engine. If the breakout time of the outflow is longer than the duration 
of the active central engine, the outflow stalls, producing no observable gamma-ray 
emission [334]. High-energy neutrinos are produced within the outflow and can escape 
through the stellar envelope, producing similar neutrino emission as for "successful" 
GRBs. Other interesting sources include low-luminosity (LL) GRBs, that can frequently 
go unobserved if their luminosity falls below the threshold of gamma-ray telescopes. LL 
GRBs, while weaker neutrino sources than their high-luminosity (HL) counterparts, 
are thought to be much more frequent, making their overall neutrino flux comparable 
or even surpass the flux from conventional HL GRBs [367, 377, 378]. Besides being 
interesting for sources not observed through their electromagnetic emission, GW-high- 
energy neutrino observations can provide information on the internal structure of the 
progenitor, as well as the properties and dynamics of the relativistic outflow (e.g., [364]). 

A potentially interesting subclass of GRBs for joint GW-high energy neutrino 
observations is low-luminosity GRBs [379, 289, 367, 380] that may be a distinct 
population from high-luminosity GRBs [289, 381, 380]. It was suggested that more 
massive stars may form a black hole after collapse, and produce a high-luminosity, high- 
Lorentz factor outflow, while less massive stars may form a neutron star after collapse, 
and drive an outflow with low luminosity and Lorentz factor [382, 381]. More recently, 
Bromberg et al. [380] found that estimated jet breakout times for low-luminosity GRBs 
are much longer than the observed durations, in marked difference with high-luminosity 
GRBs, suggesting that low-luminosity GRBs may be "choked", i.e. their relativistic 
outflow stalls before breaking out of the stellar envelope (see also [383, 384]). Low- 
luminosity GRBs are thought to be an order of magnitude more common than their high- 
luminosity counterparts, with lower electromagnetic luminosity and smaller beaming 
angle [379, 289], making them a promising source type for multimessenger detections. 

Joint searches for GWs and high-energy neutrinos are well-suited for multimessen- 
ger analyses, since both messengers have typically sub-threshold significances. For this 
purpose Baret et al. [283] developed a multimessenger analysis method that combines 
the significance and other information from the two messengers. 
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4-4- Host galaxies 

The location of potential or confirmed host galaxies of GRBs can be utilized in GW 
searches to enhance sensitivity as well as in event selection or localization. Long 
GRBs are found to occur in star-forming regions of distant galaxies, in accordance 
with their association to the deaths of massive stars [256]. Consequently, the rate of 
long GRBs is correlated to the blue luminosity of galaxies [385, 386, 256, 387, 388]. 
Long GRBs are actually found in the very brightest regions of their host galaxies, 
which are significantly fainter and more irregular than typical host galaxies of core- 
collapse supernovae, suggesting that they are associated with extremely massive stars, 
and galaxies with limited chemical evolution [387]. 

The distribution of short GRBs seem to follow that of old (~ 1 Gyr) stellar 
populations within normal star forming as well as elliptical galaxies [13, 389, 23], 
in stark contrast with the distribution of long-GRB host galaxies [13]. It is highly 
unlikely that short and long GRBs are drawn from the same underlying population 
[389]. Nevertheless, the rate of short GRBs in a given galaxy is correlated with the 
optical light of their host, and to less extent with the blue luminosity of the galaxy [390], 
indicating that short-GRB progenitors have a wide age distribution of several Gyr. As 
neutron stars may experience an initial kick at birth [391, 392, 184], binary mergers can 
take place far away from the star-forming region where they originate from. For kick 
velocities of (9(100 km s _1 ) and inspiral times of several Gyr, binaries could travel large 
distances before their merger [13, 393], traveling far outside their host galaxies. Typical 
predicted distances are ^ 10 — 100 kpc [228], the predicted distance distribution being 
well-matched by observed short-GRB distributions [394] . 

Magnetars [188, 13] form another possible progenitor type of short GRBs, making 
up less than ^one-third of their population [389] . They are weaker and softer gamma-ray 
emitters than other short GRBs, making them identifiable mostly at smaller distances, 
within the Milky Way and the Large Magellanic Cloud. Magnetars are thought to be 
created in supernova explosions [395] . As they can receive an initial kick velocity during 
the supernova explosion similarly to radio pulsars, only a fraction of them may be near 
its respective supernova remnant, in accordance with observations [395]. Magnetars are 
observed to be mostly young (^10 4 yr) objects [395], therefore they did not travel far 
from the star forming region of their respective supernovae. Consequently, extragalactic 
magnet ar flares can be expected to occur within the star- forming regions of galaxies. 

The expected distribution of GRBs can be used to enhance GW searches through 
identifying preferential GW source directions (e.g., [283, 9]). A priori information on 
source distribution has been used in multiple GW-GRB searches [226, 302, 301, 298]. 

The galaxy distribution plays a significant role in electromagnetic follow-up 
searches. Since most GW signal candidates have a reconstructed direction with an 
uncertainty of a few degrees (which may be scattered over a larger area on the sky) , the 
sky area of the possible GW source direction is greater than the field-of-view of most 
astronomical instruments used for follow-up observations. To decrease the surveyed 



CONTENTS 



35 



area on the sky, the directions in which electromagnetic follow-up is performed is down- 
selected using the directions of galaxies overlapping with GW sky area of the possible 
GW source direction [278]. 

In multimessenger searches, the distribution of galaxies can be used to increase 
the significance of astrophysical signals as well as to reject background events whose 
directions do not overlap with the direction of a galaxy (e.g., [283]). The significance 
of event candidates can be enhanced by weighing events with the a priory probability 
distribution of multimessenger sources (e.g., based on the luminosity and distance of 
galaxies) . 

Galaxy sky locations are especially important for GW sources that are detectable 
only within ~ 100 Mpc. On this distance scale, existing galaxy catalogs are more 
complete (e.g., [396]), and the number of galaxies is small within the sky area defined 
by the directional uncertainty of GW observations (~ few degrees; [9]). The use of 
galaxy direction on this scale could therefore significantly add to the sensitivity of these 
searches. Galaxy catalogs can still aid searches for sources at farther distances, although 
the sensitivity gain significantly reduces with distance [9]. 

5. Discussion & Outlook 

In this review we presented an overview of the major directions in which gravitational- 
wave (GW) astronomy has the potential of advancing our understanding of GRBs and 
their progenitors. We discussed the different promising directions that can lead to 
GW emission from GRB progenitors that is detectable with second-generation GW 
observatories. We also discussed some of the major scientific questions that could be 
addressed upon the detection of GWs from these sources. Second generation GW 
observatories, commencing their operation in the next few years (e.g., [8]), will be 
sensitive enough to reach distances for which the detection of GW signals from GRB 
progenitors becomes possible. 

The GW signal expected from compact binary coalescences, the likely progenitors 
of short GRBs, is fairly well understood as a function of the (often unknown) source 
properties. The detection of GWs can provide estimates or constraints on these source 
properties. We presented a detailed overview of the questions that can be addressed 
about short-GRB progenitors upon the detection of their GW signature. 

The core-collapse of massive stars, the likely progenitor of long GRBs, has less 
constrained theoretical GW emission models. This is partly due to the lack of our 
detailed understanding of the progenitor (e.g., how cores maintain their high rotation 
rate prior to collapse), and partly due to the complex dynamics of the collapsing core 
and the resulting compact object (protoneutron star or black hole) and accretion disk, 
whose understanding requires further theoretical and numerical studies. Nevertheless, 
the core-collapse dynamics may lead to GW signals that are detectable on distance 
scales comparable to the lower end of the long-GRB distance scale (> 100 Mpc). For 
example rotational instabilities in differentially rotating protoneutron-stars, as well as 
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global instabilities in tori around black holes, both lead to rotating, non-axisymmetric 
systems with potentially strong GW emission. 

Starquakes in highly-magnetized neutron stars, likely responsible for producing a 
fraction of short GRBs, result in seismic oscillations in the neutron star, which in turn 
may produce GWs detectable with second generation GW observatories. The detection 
(or non-detection) of such GWs coincident with the flaring activity of neutron stars can 
be informative, e.g., on the properties of matter at nuclear densities, as well as on the 
properties of astrophysical neutron stars. 
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